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Abstract 

Several results on the reconfigurable control architec- 
ture for the formation flying of multiple spacecraft 
are presented. In this direction, simple control laws 
are combined with logic-based switching to propose 
a hybrid control architecture for leader reassignment, 
leader-following capturing, and dealing with control 
saturations. 
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1 Introduction 

Formation flying (FF) has been identified as an en- 
abling technology for many of the NASA’s 21 st cen- 
tury missions, among them, the Space Technology- 
3 (ST-3) and the Terrestrial Planet Finder (TPF). 
Formation flying involves flying a group of space- 
craft in a particular pattern while maintaining precise 
(possibly time varying) relative position, velocity, at- 
titude, and angular velocity, with respect to each 
other [3], [6]. Since traditional spacecraft control is 
often concerned with measuring and maintaining the 
same quantities for a single spacecraft with respect 
to an inertial reference frame, the analogous FF con- 
trol and estimation problems are often an order of 
magnitude more challenging than those encountered 
traditionally for a single spacecraft [1], [4], [8]. In 
order to make the FF control problems at least sim- 
ilor to the single spacecraft case, an approach based 
on leader- following has been proposed by Wang and 
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Hadaegh [7]. The basic idea in leader- following is 
to designate a particular frame (or multiple frames) 
in the formation as the reference frame(s) of inter- 
est and measure and control the states of the rest 
of the formation with respect to them. The purpose 
of the present paper is to show that linear matrix 
inequalities (LMIs) [2] can be combined with logic- 
based switching schemes to propose reconfigurable 
control architecture for the formation flying of mul- 
tiple spacecraft. 

The outline of the paper is as follows. In §2 the no- 
tation used in the paper is presented. Simple control 
laws for the formation flying control are then derived 
in §3 based on the leader- following concept. In §4, 
§5, and §6, the control laws derived in §3 are com- 
bined with logic-based switching to propose a hybrid 
control architecture for leader reassignment, leader- 
following capturing (defined subsequently), and deal- 
ing with control saturations. 

2 Notation 

Formation flying consists of flying a group of space- 
craft in a particular pattern. To be able to express 
the time evolution of the formation and design the 
corresponding control laws, it is convenient that a 
reference frame is attached to each spacecraft. We 
shall always assume that these reference frames are 
induced from a dextral set of three orthonormal vec- 
tors. Let the formation have n spacecraft labeled as 
l,2,..,,n. Let T x denote the reference frame at- 
tached to the z-th spacecraft; j F* on the other hand 
shall designate the inertial reference frame. For the 
inertia and the mass of the i-th spacecraft we use /* 
and m x , respectively. The force and torque acting 



Figure 1: Formation Coordinates 


3 Simple Control Laws for Leader- 
following 

In this section we go over some simple control laws 
for formation flying that are derived based on the 
state feedback synthesis procedure which uses LMls 
as its building block. These control laws can be used 
for the control of the formation pattern under two 
different measurement scenarios. First, we consider 
the situation where inertial measurements are avail- 
able to both the leader(s) and the follower(s); then 
we comment on the case where the follower(s) mea- 
surements are done with respect to its own moving 
reference frame. In all subsequent sections, we say 
that i is the leader of j if is an affine function of 
r 1 which is twice differentiable. 


upon i are denoted by /* and r* ; for the mass nor- 
malized force we used u‘ := -4-. The time derivative 
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with respect to shall be denoted by likewise, 

37 be used for the time derivative with respect 
to T 1 . denotes the position of the origin of P 
with respect to P ; r l is the position of the origin of 
P with respect to P . The desired position of the 
origin of P with respect to P shall be denoted by 
T <i > by r d when j = I. The velocity of the origin 
of P with respect to P , the velocity of the origin 
of P with respect to P , the desired velocity of the 
origin of P with respect to P , and the desired ve- 
locity of the origin of P with respect to P , shall be 
denoted by , v l , , and v' d , respectively. The 

vector [r* v l ] r shall be referred to as the state of the 
i-th spacecraft and will be denoted by x 1 . Similar 
notations are used for the attitude and the angular 
velocity of P with respect to P : and are the 

attitude and the angular velocity of P with respect 
to P and q d and u d are the desired angular veloc- 
ity and attitude of P with respect to P (refer to 
Figure 1). All other notations are standard: 9i n de- 
notes the real Euclidean space of dimension n; H-Hoo 
and ||.|| are used for the infinity norm and the 2-norm 
for vectors and matrices. The cross product matrix 
induced by the vector x = [x x x 2 x 3 ] T is the matrix, 


3.1 Inertial Reference Frame Measure- 
ments 


Let i be the leader of j during the time interval [t 0 , f/]. 
The desired position of j is thus expressed as, 

rj(t) = r*(f) + h u (t), t 0 <t <t f . 


The error expression for j is then simply, 

= r d(t) ~ ^(0 = r j (0 - H(t) + /»«(*). 


Assuming that h" is twice differentiable on [*„,*/], 
the above expression can be differentiated twice with 
respect to the inertial reference frame to obtain (re- 
calling that u* = P), 


dt 2 


= u l (t) - u>(t) + 


d 2 hV(t) 

dt 2 ' 


By letting, 
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= u* + 


d 2 hv(t) 

dt 2 


+ **(0. 


one obtains. 


d 2 e>(t) 

dt 2 


= -«*(«) 


(3.2) 
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(i.e., feedback linearization). The equation (3.3) can 
be expressed as, 
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where *i(f) _ e^(f), z 2 (t) = , and the matrices 

/l J and B ^ are defined as suggested by (3.4). 

The control design based on the state feedback syn- 
thesis for leader- following is thus reduced to finding 
the term for u j using the following LMI, 

AiQ + QiA *) 7 + B*Y + Y T (B i f < 0 , ( 3 . 5 ) 

Q > 0 , ( 3 . 6 ) 

and let K* = YQ~ X [2]. Hence, given that i is the 
leader of j the control law for j has the form, 

«*(*) = «‘(0 + - + rg-V(0, t 0 <t< 

Employing the above control law by the follower 
spacecraft j guarantees that the origin is the globally 
asymptotically stable equilibrium of the error func- 
tion z(t), and thereby, r^(f) -+ r^(f) as f -*■ oo. 

3.2 Moving Reference Frame Measure- 
ments 


*a(0 = -2z 3 (t)x z 2 (t)-(I i )-\ T i(t) 

~ * 3(0 X f‘* 3 ( 0 ) - 23(0 x (* 3 (f) x * x (f)) 

+ (u*(f) - vf(t)) + - ^ 1 , ( 3 . 10 ) 

* 3(0 = (/ i )~ 1 ( rJ (0 ~ 23(0 x Pz 3 (t)). (3.11) 

Consider two distinct situations. 

1 - j has constant angular velocity: Consider the 
case where, 

r j (t) = * 3 (f) x Pz 3 (t), ( 3 . 12 ) 

i.e., the angular velocity of j during the leader- 
following remains constant. The dynamical 
equations ( 3 . 9 )-( 3 . 10 ) can then be written as, 

2l(0 = 2 3 (f), 

22(f) = Wizi(t) + W 2 z 2 (t) + u ! (f) - u-i(f), 
where, 


We shall now briefly go over the situation where the 
measurements are done in the moving frame attached 
to the follower spacecraft. Feedback linearization is 
then used to reduce this case to that considered in 
the previous section. 

Again let i be the leader of j during the time in- 
terval [f 0 ,f/]. Contrary to the case considered pre- 
viously, we would like to obtain an expression which 
describes the time evolution of e 3 in T 3 (as opposed 
to F l ). Proceeding from (3.1), one has, 


W(f) x (J(f) x e*(f)) 

= (u ! (f) - u 3 (t)) + - (3.7) 


the last term on the right hand side of (3.7) can of 
course be represented in T 3 . 

The rate of change of the angular velocity uP with 
respect to J- ' 3 (or jA ) is related to the applied torque 
on the spacecraft via the Euler’s equation, 


dJ{t) 

dt 


= - uA(t) x 


( 3 . 8 ) 


Let *,(f) = e>(<), z 2 (t) = and * 3 (f) = ^(f). 

The dynamics of j can thus be represented as, 


2t(<) = z 2 (t), 


^1 = 232^-11*3112/, and W 2 = -[2* 3 ], 

Consider again the change of variable of the form 
u 3 {i) = u‘(f) + -(- «j(f); then, 
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Define the matrices A 3 and B 3 as suggested 
above; we can now proceed to solve the LMI, 

+ q{a 3 ) t + RJy + y t {b 3 ) t < o, (3.13) 

Q > 0, ( 3 . 14 ) 


and let, 

”*(0 = u*(f) + -- + YQ-'P(t), f 0 <f<f/; 

note that only the definition of the matrix A* 
has been modified from that used previously to 
reflect the fact that the error vector is now mea- 
sured in the moving coordinate frame attached 
to the follower. 

2 . j has non-constant angular velocity: If the angu- 
lar velocity of j does not remain constant during 
the leader- following, then we can use feedback 
linearization to linearize the dynamics in such a 


( 3 . 9 ) 



way that the LMI approach above can still be 
adopted. For this purpose it suffices to let, 

a j = -2 z 3 (t) x z 2 (t) - (/ i )“ 1 (rJ( t) 

-z 3 (t) x /*z 3 (0) ~ *3(«)(*3(<) x *i(0), (3.15) 

and let u j (t) = u^f) + — + u j (f) + aty); as 
before the expression for £•*(/) is found be solving 
an LMI. 


In both scenarios considered above, the control law 
for the leader spacecraft i can also be based on the 
state feedback synthesis. For this purpose it suffices 

to let u l (t) = where the matri- 

ces Y and Q are found from the LMI (3.5)-(3.6) by 
letting, 


A l = 


0 I 

0 0 ’ 


and B l = 



however z is now simply rfa) - r‘(t). 


4 Leadership Re- Assignment 

The designation of the leader, aside from its asso- 
ciated hardware and software considerations and the 
required communication protocol, is rather arbitrary. 
It is thus of interest to consider a situation where the 
leader assignments are time varying, and that any 
subset of the spacecraft in the formation can assume 
the leadership role. 1 

Suppose that at a particular instance of time, i 
is the leader of j; in this case the control law of §4 
(when inertial measurements are available) can be 
implemented as, 

u l (t) = Kz\t)+^l, 

**(0 = A'ri(t) + u'(t) + 

where z l (t) is the state error observed by i at time 
t. Since h u (f) = — hJ*(f), when the leadership assign- 
ment is reversed and j is the leader of i, the control 
laws can be reconfigured as, 

= AV(0 + »l(0 - 

«‘(0 = AV (!) + ^5^; 

1 Leadership reassignment becomes specially relevant when 
one looks beyond ST-3 and TPF type space interferometry 
missions to the formation flying of large number of spacecraft. 
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Figure 2: Switching for Leader Reassignment 
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Figure 3: Leader- following capturing 


refer to Figure 2 

5 Leader-Following Capturing 

We consider a situation where a free spacecraft is 
captured by a leader following scheme. Building on 
the control laws developed in §4, the corresponding 
block diagram representing the switching control sys- 
tem can be drawn as in Figure 3. Note that we have 
considered the situation where the isolated spacecraft 
is not assigned as a leader; if this is in fact the case, 
then its control law will not be changed from when 
the spacecraft was free. However, in this latter situa- 
tion, the control law for the new followers of the new 
leader spacecraft changes according to the procedure 
presented in §5. 

6 Control Saturations 

We now examine the scenario where the j-th space- 
craft, j, following i in a leader-following scheme, is 







also avoiding control saturation by switching between 
two or more controllers. In order to simplify the pre- 
sentation, we shall assume in the rest of this sec- 
tion that — =0, to < t < tj. Let m denote 
the 2-norm of the maximum allowable mass normal- 
ized force on each spacecraft in the formation; that 
is, we require that, ||t*l(t)|| < m, t 0 < t < t f . 
Now since iP(t) = YQ~ 1 z(t ) + u‘(t) we require that, 
lk(0 + YQ~ l z(t) || < m, t 0 < t < t f . Note that 
although j has no prior knowledge about the values 
of u‘ it has to choose Q and Y such that the con- 
trol constraint is satisfied. To cope with this lack 
of knowledge on the values of u l (t), we proceed to 
present a controller switching mechanism which sat- 
isfies the control constraint, in face of the lack of 
a priori knowledge of the values of u\t) by the fol- 
lower spacecraft. The only assumption which is re- 
quired for the proposed approach to work is that, 
IK(0ll < m i to < t < tf. Let us start with the 
stronger requirement, \\YQ~ l z{t)\\ < m - ||u i (/)|| = 
m* (t), to < t < tf in order to satisfy the control 
constraint. Let S to = {z : z T Q; o 1 z < 1}, where Q h is 
a positive definite matrix which is chosen such that 
^(0) belongs to £ to by solving an LMI, in conjunction 
with, 


AQt 0 +Qt 0 A+ BY t0 +Y*B t <0. (6.16) 

For small values of St, if x{(to) = x{(to + St) for t <E 
[*o,*o + &t] and we use the controller K to = Y tg Qt 0 l , 
then it would be the case that z{Iq + St) 6 S to . In 
fact, if arj,(t) remains constant, then z(t) € S to for all 
t € [to,*/]- In this situation, in order to guarantee 
that the saturation constraint is not violated, we can 
augment the above LMIs with another one, 



Figure 4: Ellipsoids for Control Switching 


AQt 0 + Qt 0 A T + BY t0 + Y t T 0 B T < al, 
Qt 0 >0, a < 0, 

f i * r (o) 

I *(0) Q t0 


Qto YZ 

Y t0 m‘ 


> 0 , 

> 0 , 


We then proceed from time t 0 and considered the 
various scenarios which can occur at time to + St: 

1. € £t 0 and m l (t) has remained constant. 

2. z(to + St) € £ to 1 however m l has changed over 
the interval [/ 0 , *o + tf/]. 

3. z(t 0 + St) £ £ to , whether or not m 1 has remained 
constant. 


For each scenario above, we were able to show that 
an LMI can be solved to address the control satura- 
tion problem [5], Moreover, the resulting proposed 
switching mechanism results in a hybrid dynamical 
system for which the origin is the globally asymptot- 
ically stable equilibrium point. 


7 Conclusion 


' Qt 0 Yg 

Y to m\t) 


> 0, 


since [2] max<> 0 ||u j (t)|| = maxt> 0 ||F (o Q^ 1 2 (t)|| < 

max iefto ||y t0 Q- I .(«)|| < \ m t X (Q- l/ 2 Y t T o Y to Q- l/2 ). 
The problem is that in general, one cannot guar- 
antee that z(to + St) 6 £ to , nor does the above 
discussion addresses the situation where m‘(£) does 
not remain constant. We are thus led to incorpo- 
rate logic-based switching in conjunction with above 
LMIs to address both of these scenarios. Let m* := 
ll ml (0H; solve the semi-definite program, 


min 

r t0 


Several reconfigurable control strategies for the for- 
mation flying of multiple spacecraft were presented. 
In this direction, it was demonstrated that by em- 
ploying feedback linearization and linear matrix 
inequalities, in conjunction with simple switching 
mechanisms, various formation control and manage- 
ment issues could be addressed. 
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